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Abstract 

We compute the thermal conductivity and sound attenuation length of a dilute atomic Fermi gas in the 
framework of kinetic theory. Above the critical temperature for superfluidity, T c , the quasi-particles are 
fermions, whereas below T c , the dominant excitations are phonons. We calculate the thermal conductivity 
in both cases. We find that at unitarity the thermal conductivity k in the normal phase scales as k oc T 3 / 2 . 
In the superfluid phase we find k k T 2 , At high temperature the Prandtl number, the ratio of the 
momentum and thermal diffusion constants, is 2/3. The ratio increases as the temperature is lowered. 
As a consequence we expect sound attenuation in the normal phase just above T c to be dominated by 
shear viscosity. We comment on the possibility of extracting the shear viscosity of the dilute Fermi gas at 
unitarity using measurements of the sound absorption length. 



I. INTRODUCTION 



Cold, dilute Fermi gases in which the interaction between the atoms can be tuned using an exter- 
na> magnetic field provide a aew paradigm for strongly correlated quantum fluids HQ. Recent.y, 
there has been significant interest in the transport properties of these systems [3J. This interest was 
triggered by three independent developments. The first was the experimental observation of almost 
ideal hydrodynamic flow in the normal phase of a cold Fermi gas in the limit of infinite scattering 
length 4|. The second was the discovery of almost ideal flow in a completely different system, the 
quark gluon plasma created in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) 
5KZJ. Both of these experiments constrain the shear viscosity of the fluid. Remarkably, while the 
absolute value of the shear viscosity of the two systems differs by about 26 orders of magnitude, 
the ratio of shear viscosity rj to entropy density s is very similar, rj/s ~ (0.1 — 0.5) h/ks, where h 



is Planck's constant and ks is Boltzmann's constant [8l-ll2|. 



The third development was the theoretical discovery of a new method to compute the shear 
viscosity of strongly coupled field theories [3, Q. This method is based on the hologra phi c 
duality between four dimensional field theory and string theory in higher dimensional spaces 15]. 
Using the holographic equivalence one can show that in a large class of field theories the strong 
coupling limit of t]/s is equal to h/(4irkB) — 0.08 h/ks- This value is remarkably close to the 
experimental results for rj/s quoted above. Initially the value h/^^ks) was derived only for scale 
invariant relativistic field theories, but it was later shown that the same value of rj/s is obtained in 
the strong coupling limit of many non-conformal and non- relativistic field theories 16l. Il7|. 

In this work we will use kinetic theory to compute another important transport property of a 
dilute atomic Fermi gas, the thermal conductivity k. We will consider both the high temperature 
as well as the low temperature, superfluid, phase. Our results complement earlier studies of the 
shear and bulk viscosity in the high and low temperature phase. The shear viscosity was studied 
in In the unitarity limit the bulk viscosity of the high temperature phase vanishes [2oj|. In 

the superfluid phase there are three bulk viscosities, one of which is non-vanis hing in the unitarity 
limit. The bulk viscosities in the superfluid phase were recently computed in [2l|. Combining our 
results for the thermal conductivity with the existing calculations of the shear viscosity we can 
study the relative size of momentum diffusion and thermal diffusion. The ratio of these quantities 
is known as the Prandtl number 

Pr=^, (1) 
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where cp is the specific heat at constant pressure and p is the fluid mass density. 

The Prandtl number also controls the relative importance of shear viscosity and thermal conduc- 
tivity in sound attenuation. Measurements of sound attenuation have the potential to provide new 
experimental constraints on transport properties of dilute atomic gases. Current estimates of the 
shear viscosity are based on the study of scaling flows, such as collective modes or the expansion out 
of a deformed trap [10;]. These experiments do not directly constrain the density independent part 
of the shear viscosity. This is a problem, because kinetic theory predicts that the shear viscosity 
at high temperature is only a function of temperature, and not of density. Sound attenuation, on 
the other hand, is sensitive to this contribution. Using our results for the Prandtl number it is 
possible to use measurements of sound attenuation to constrain the shear viscosity. Conversely, in 
the regime where shear viscosity is well constrained it is possible to extract thermal conductivity 
from sound attenuation. First sound was observed in j^ . and ideas for detecting second sound are 
discussed in [23l |. 

This paper is organized as follows. In Section UU we will review the quasi-particle interactions 
in dilute atomic Fermi gases, and in Section IIHI we outline the calculation of thermal conductivity 
using kinetic theory. In Sections IIVI and IVl we compute the Prandtl number and sound attenuation 
length. We end with a few comments on the possibility of extracting the sound attenuation length 
from experiments with optical traps in Section IVII Appendix [A] contains supplementary material 
regarding the trial functions used for solving the linearized Boltzmann equation, and in Appendix 
[B]we collect some thermodynamic relations. 

II. ELEMENTARY EXCITATIONS 

Kinetic theory is based on the existence of well defined quasi-particle excitations. In the dilute 
Fermi gas quasi-particles exist both at very high and at very low temperature. In the high temper- 
ature limit the relevant degrees of freedom are spin 1/2 non-relativistic atoms with mass m. In the 
dilute limit the scattering cross section is dominated by elastic s-wave scattering. The scattering 
amplitude can be described in terms of a delta-function interaction V(x — x') = C 5(x — x 1 ), see 
Fig. [TJ Here, C = Ana/m is related to the scattering length. The s-wave cross section is determined 
by sum of all two-body diagrams. The cross section is 




(2) 
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FIG. 1: This figure shows the main processes that determine the thermal conductivity in high and low 
temperature phase, respectively. Fig. (a) shows the local four-fermion interaction, and Fig. (b) shows the 
three-phonon interaction. 



where Pi,P2 and P3,P4 are the ingoing and outgoing momenta and q is the relative momentum 
q = (pi — p2)/2. The relative momentum of the outgoing particles is q' = (p3 — P4)/2 with 
|q'| = |q|. The solid angle is defined in terms of the angle between the relative momenta q and q', 
dQ = G?cos(# qq ')<i0 qq /. 

The characteristic temperature of the system is the Fermi temperature Tp = hp/ (2m), where 
kp is the Fermi momentum defined by kp = 37T 2 n to t and n tot is the density of both spin states. 
The strength of the interaction can be characterized by the dimensionless parameter (kpa). The 
regime of small and negative (k F a) is the BCS limit. In this limit the Fermi gas becomes a BCS 
superfluid at a critical temperature T c which is exponentially small compared to Tp. The ratio 
TjjTp increases with |/cpa|. In the unitary limit \kpa\ — > 00 the critical temperature is T c ~ 0.15Tp 



241 ]. The regime of positive (hp a) is called the BEC limit. For small and positive (kpa) the system 
is composed of weakly interacting bosons with mass 2m and density n tot /2. Superfluidity occurs at 
the Einstein temperature T c ~ 0.21Tp. 

In the superfluid phase the low energy excitations are phonons, Goldstone modes associated 
with the spontaneous breaking of the U(l) phase symmetry. At unitarity the effective lagrangian 



for the phonon field is strongly constrained by scale invariance. Son and Wingate showed that 



25|, 



= \W) 2 - ^ S 2 (V0) 2 - g 3 [(<9 O 0) 3 - 9t; s 2 «9o0(V0) 2 ] (3) 
-\gl [(<9o0) 4 + 18*; 2 (<9 O 0) 2 (V0) 2 - 27v' s ( V0) 4 ] + . . . , 

and v s is the speed of sound given 



where g 3 = 7rt' 3/ ' 2 ^ 3//4 /(3 1 / 4 8 / u 2 ), /i is the chemical potentia 
by v 2 s = 2/x/(3m). At unitarity \i = ^E F with ^ ~ 0.44 
dispersion relation is linear, e p = v s p. We will see that the thermal conductivity is sensitive to 



261 ] . At low momentum the phonon 



non-linearities in the dispersion relation. We will write the dispersion relation as e p = v s p(l + 7P 2 ), 
where the parameter 7 can be expressed in terms of the coefficients of certain higher order terms in 
the lagrangian of Son and Wingate, 7 = — 7r 2 \/2£(ci + l.hc-zj/k'p. Rupak and Schafer computed c\ 
and C2 using an epsilon expansion [27[. They find c\ ~ — \/2l (307r 2 £ 3//2 ) and c 2 ~ 0, corresponding 



to 7 = l/(30m/i). Other estimates of 7 can be found in 28430 1. 

We will study the contribution to thermal conductivity from the phonon splitting process ip ^ 
(p + (p, see Fig. [TJ For 7 > this process is kinematically allowed. The matrix element for the three 
phonon process is 

iMz = g 3 (p [Q ■ K] + k [Q ■ P] + q [P ■ K]) , (4) 

where we have defined the four vector product Q ■ K = q^k^ — 9v^q ■ k. If 7 is negative then the 
dominant contribution arises from (p + (p <p + (p. We recently studied this process in a relativistic 



superfluid, the color-flavor- locked phase of QCD [3JJ. We found that the 2 ^ 2 processes are 
controlled by s, t, w-channel exchanges with almost on-shell intermediate state phonons. This means 
that the dominant process can be viewed as 1 2 phonon splitting followed by 2 ^=f 1 phonon 
absorption. 

III. THERMAL CONDUCTIVITY 

Thermal conductivity is defined by the relation between the heat current h and the temperature 
gradient, h = -kVT. In kinetic theory the heat current is expressed in terms of the quasi-particle 
distribution function, 

h = v J dTe p v p 5f p (5) 

where v is the degeneracy factor, dr = d 3 p/(27r) 3 , v p is the particle velocity and 5f is the deviation 
of the distribution function from its equilibrium value, 5f = f—f°, with /o(p) = l/(e^ p_At ^ T ±l) for 
bosons and fermions, respectively. Note that in the case of phonons there is no conserved particle 
number and \x = 0. The correction to the distribution function is linear in VT and we can define 

With this form we can rewrite Eq. (JHJ) as 

z/VT 



h 



3T 2 



J dTe p v p ■ x/o(p)(l ± kip)) = -TVT(e p v p \ X ) , (7) 
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where we have defined an inner product as 

(a\b) = ^5 JdTf (p)(l±f (p))a-b. (8) 
We can then write the thermal conductivity as 

K = T(e p v p \ X ). (9) 

There are certain constraints on the the form of X- The Boltzmann equation conserves energy, 
momentum, and, in the case of a conserved charge, particle number. A non-trivial constraint on x 
arises due to momentum conservation. We have 

= vj dT P 5f p = -^J dl>(x • VT)/ (p)(l ± / (p)) = -TVT(p\x) , (10) 

which implies that (p\x) = 0. As we can see from the constraint and Eq. ([9]), any term in e p v p 
that is linear in p will not contribute to the thermal conductivity. In particular, if the kinetic 
description is governed by a quasi-particle with exactly linear dispersion relation then the thermal 
conductivity vanishes. 

The function x is determined by solving the Boltzmann equation 

| = f + !V V/ = C[/], (11) 

where C[f] is the collision integral. Schematically we can write 

C[f) = J dr ft _iiu(m; out)D nMout , (12) 

where dr n _i is the (n — 1) particle phase space (all the particles not labeled by momentum p), w is 
the transition probability, and D in ^ ont contains the distribution functions. Taking the convective 
derivative on the left-hand side of the Boltzmann equation and focusing only on terms proportional 
to the temperature gradient we can write 

d/ / (l±/o) _ /o(l ± f0) I \ ^jrr, /i q\ 

■ = — a v ■ Vi = — \a v ) ■ vi (IS) 

dt 3T p 3T 1 vl K J 

where cx p will depend on equilibrium properties of the quasi-particle system. In the case of either 
particles with linear dispersion (like phonons) or particles with non-relativistic dispersions and a 



chemical potential, we can write, [3JJ , 



^ ~ ^ + Ap (14) 



where A depends only on the thermodynamics of the system. We can linearize the right-hand side 
of Eq. (ITTI) in x as 



C[f] = C[f ] + C[6f] = C[Sf] = J dVw(in; out)5D irMout = /o(1 ^ /o) C\ X ) • VT , (15) 

where we have defined the linearized collision operator, C, acting on the state \x)- We will specify 
C in more detail below. We can now write the linearized Boltzmann equation as 

\cx p ) = C\ X ). (16) 

The thermal conductivity is now determined by solving Eq. (1161) for \x) and then computing k 
using Eq. (Q. In practice, we will adopt a variational procedure. Using the constraint we can show 
that T(x\ct p ) = (6pV p \x) for both fermions and bosons. The thermal conductivity can then be 
written in two alternative ways 

K = T(e p v p \ X ) and K = T a {x\C\ X ), (17) 

where we have used the Boltzmann equation, Eq. (IIEI) . to derive the second equality. These two 
representations form the basis of the variational principle. Consider a trial state \g). The Schwarz 
inequality implies that 

(g\c\g}(x\c\x}>\(9\c\x}\ 2 . (is) 

Using Eq. (ITT)) and Eq. (fl6)h we can write this as 

iz>T 2 \(9\<*p)\ 2 _ \(9\e P v P )\ 2 (1Q) 
~ (g\C\g) (g\C\g) 

Using Eq. (ITT)) we observe that an exact solution of the linearized Boltzmann equation saturates 
the inequality. We can then choose to expand the trial function g in orthogonal polynomials and 
take advantage of that fact that g is a dimensionless vector-valued function to write 

N 

g = X J2bi F > H) Bi F > H \x>) (20) 

where x = p/ A, x = \x\, and A is some characteristic energy scale in the problem (e.g. A = T/v s 
for massless particles and A = y/2mT for massive particles). The polynomials are such that 
B = l,Bi = x 2 + a 10 , B 2 = x 4 + a 2 ix 2 + a 2 o, etc., and (F,H) labels the basis functions in the 
fermion and phonon case. The orthogonality condition is 

Ai F ^S st = J dxf (x)(l ± f (x))x*Bi F ' H \x 2 )B { t F ' H \x 2 ) . (21) 



The type of quasi-particle involved enters via the statistical factor (1 ± fo(x)) and through the 
dispersion relation in fo(x). The specific forms of Bm' H \x 2 ) are given in Appendix |A] Note that 
we have enforced the momentum constraint by excluding the s = term in Eq. (j20p . 



A. Normal Phase 



In this section we will specify the matrix elements and the collision operator in the case of 
atom-atom scattering in the normal phase of the dilute Fermi gas. We will scale all the momenta 
aS X{ = Pi use the Fermi-Dirac distribution functions and set v = 2. The numerator of 

Eq. (TT9D is 

(9\e P v p ) = 3^3 / dr/ o (p)(l-/o(p))e p tV0 (22) 
= S^XX / dx/ (x)(l-/o(x))x 6 5f(x 2 ) 
Am sr^ h F A F, _ Am h F a f 

where we have used the orthogonality condition. The matrix element of the collision operator is 

(a\C\g) = ^ X W / dn l«i " ^1 YVaW D ° 2 ^ 91 ' [Ql + 92 ~ 9s - 9 * ] ■ (23) 

Here we have written the transition probability in terms of the cross section Eq. (j2]) and the flux 
factor \vi — v 2 \. The integration variables are the pair momentum P = p Y + p 2 , the relative 
momentum q = (p 1 — p 2 )/2, and the solid angle Q is defined below Eq. (j2j). We have also defined 
D 2^2 = fo(x 1 )f (x 2 )[l - fo(x 3 )][l - foM]. We can write 

( n \r\n\ 2 n / 2 m 3 / 2 ^ F F f d 3 xd 3 y mTa 2 y F F 

(g\C\g) - -^-ij.bt J -j^yr^ i + 2mTa2y2 D 2 ^ 2 B s B t (24) 

911/2,^3/2 
_ I'm' b F b F M F 

~ 32^3/2 ° s °* lvlst ' 

wherex= \P\/y/2mT,y = \q\/V2mT and B( = [B[(x 1 )x 1 +B[(x 2 )x 2 -B[(x 3 )x 3 -B[(x 4 )x A ]/2 
. We finally get 

3 

kf > {9h^f = (T_y Dp {kpa) t/Tf) j (25) 



(9\C\g) 3 • 2V%4 VTf 
where we have defined and 



D ^' TIT ^Smk- (26) 
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FIG. 2: Thermal conductivity in the normal phase of a dilute Fermi gas in units of kq = m ' 2 Tp /2 . In the 
left panel we show k/kq as a function of T/Tp for different values of kpa. In the right panel we show k/kq 
as a function of kpa for different values of T/Tp. 

We can now consider a finite basis and optimize the inequality with respect to the coefficients bf . 
We obtain Dp (kpa, T/Tp) = (A F 1 ) 2 (M F )^ where (M F )^ refers to the (1,1) -component of the 
matrix inverse of M. Note that this result is equivalent to using the matrix form of the Boltzmann 
equation in a finite basis. The function Dp must typically be calculated numerically. We can 
analytically evaluate Dp in the limit of high temperature T/Tp 3> 1 and by keeping only the 
leading term in the polynomial expansion. We get 



I \l m a 2 \®\ ^ ^ /c\n\ 

k f = < v (27) 

1 3 

=maT2 I a I — » oo 



128^ 

We have tried to improve this result using higher order terms in the polynomial expansion, but 
the correction is very small, 5k/ k < 2%. This is consistent with the results in }32| . where the 
authors observed very rapid convergence for the shear viscosity in the high temperature phase. 
Numerical results are shown in Fig. [2j Note that the thermal conductivity depends only on the 
square of the scattering length so we plot the thermal conductivity as a function of fcp|a|. We 
observe that, qualitatively, the temperature dependence of k is the same for all values oi kp\a\. In 
the unitarity limit k, grows as ~ T 3//2 for T 3> Tp. In the BCS limit the temperature dependence is 
k ~ T 1 / 2 . The thermal conductivity has a minimum at T/Tp ~ (0.2 — 0.3), and increases for very 
small temperatures. This rise is related to Pauli blocking. We should note, however, that except in 
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the BCS limit the results in this section are not reliable for T <C Tp. The right panel shows that 
at fixed temperature k decreases as the strength of the interaction increases, and that k becomes 
insensitive to kpa as we approach the unitarity limit kpa — > oo. 



B. Superfluid Phase 

In this section we study the thermal conductivity in the superfluid phase. We will concentrate 
on the unitarity limit kpa — > oo. In the BCS limit the critical temperature is exponentially small, 
and the thermal conductivity is dominated by fermionic quasi-particles even in the case T < Tp. 
n the extreme BEC limit the thermal conductivity reduces to that of an almost ideal Bose gas 



331 ] . Following the procedure outlined in the case the normal phase, we can write the numerator 



of Eq. (TI9D as 

(9\e P v p ) = ^ j dVf (p)(l + f (p))e p v p -g (28) 



2 " 7 ' 3 J>f / dxf (x)(l • f (x))x«B?(x 2 ) 



3tt 2 v* 



27 



H A H 



where we have used v = 1, scaled the momenta as x = vp/T and kept only the leading order terms 
in the parameter 7 that governs non-linearities in the dispersion relation. We note that the result 
is indeed proportional to 7, and k vanishes in the case of an exactly linear dispersion relation as 
discussed before. The matrix element of the collision operator is 



(9\C\g) = ^4 E W/ dT pkq \M 3 (p,k, q )f(27rmP + K-Q)DUB» ■ B? (29) 
£&ffc? J dT xyz \M 3 (x,y,z)f (2n)W(X + Y - Z) DUB" ■ B? 



96tt3 r<; 
^3 1 i.HuHnsrH 



E W M i 



96tt 3 v e s " s " st ' 

where dT pkq is the 3-particle Lorentz-invariant phase space, X = (x ,x) = (p /T,v s p/T),etc, 
Dl« 2 = f (x)f (y)(l + /„(*)) and Bf = [Bf (x)x + B?(y)y - Bf{z)z]/^. As before, we can 
define 

D B = {A*)*(M*)£ , (30) 
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FIG. 3: Variational function Dh as a function of the number of terms in the polynomial expansion of 
g(x). We show Dh for two different temperatures. 



which allows us to write the thermal conductivity as 



. (9\e P v P ) 2 128 7 2 T 2 
kh > , , a,, . = t; 9 n H 



25tt 3 £ 2 °\T F J V ' ' 



(31) 



(g\C\g) 3tt gi vl 

where £ ~ 0.44 is defined in Section [Til Dh is a dimensionless function of T/Tp which has to be 
calculated numerically. The convergence of Dh with the number of terms in the trial function is 
shown in Fig. [3j We observe that the result in converging, but not nearly as fast as it did in the 
case of fermionic quasi-particles. The temperature dependence of Dh is very weak and kh ~ T 2 . 
In Fig. H] we combine our results for the thermal conductivity of the dilute Fermi gas at unitarity in 
the high and low temperature limits. The dashed line indicates the location of the phase transition. 
We note that the low and high temperature results do not match very well near T c . This suggests 
that additional processes or excitations must become relevant in this regime. We also note that at 
T c critical fluctuations are expected to lead to a divergence in the thermal conductivity 34]. 



IV. SHEAR VISCOSITY AND THE PRANDTL NUMBER 



A. Shear Viscosity 



In order to compare the magnitude of thermal and momentum diffusion in dilute atomics gases 
we also need to determine the viscosity of the gas. The shear viscosity of the dilute Fermi gas at 
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FIG. 4: Thermal conductivity k/kq of the dilute Fermi gas at unitarity as a function of T/Tp. We show 
the calculations in the high and low temperature limit, extrapolated to T ~ T c . The position of the critical 
temperature is indicated by the dashed line. 



unitarity has been calculated in both the high temperature [18] and low temperature limits [19 ]. 
Here, we quickly summarize the results and present a simple numerical formula that can be used 
to compute the contribution from binary scattering of atoms for all values of k F a and T/T F . 
The viscous correction to the stress tensor is 

8Tit = uj dTv p ^8f p = -rjVii , (32) 

where Vij = diVj + djVi — 2/35^' V ■ V and Vi is the velocity of the fluid. We use an ansatz for Sf p 
analogous to the one that appeared in the calculation of the thermal conductivity, 

Sfr= _hm±iMi XijV%] . pa. 

The shear viscosity 77 is then given by 

V = J dT/o(p)(l ± fo{p))^PiPjXij = —(PijlXij) , (34) 

where pij = v p (piPj — ^SijP 2 )/p and we have used the inner product defined previously. The 
linearized Boltzmann equation for a fluid perturbed by a small shear stress is given by 

d/ = /.&,)(! ± /.&,)), } 

dt 2T 1 3/ 3 y ' 

= c[h ^sjmpMic\ Xil)Vil . 

11 



As before, we can combine these results into a form that is suitable for a variational treatment. We 
write 

(XiJ -b~,r 

5 <X«|C|x«> 

and Xjj = ^g(x). We expand g(x) = Yl s =o c s C s (x 2 ) where the polynomials C s (x 2 ) are defined in 
analogy with Eq. (l2Tj) but with the weight factor x 4 replaced by x 6 , see Appendix |A] In the case 
of binary atomic scattering the variational bound is 

o7/2 / j> \ 3/2 

>7f > -r^ 7^ ^(foro, T/Tp) , (37) 



45 VTiV 

where ??o = (mT F ) 3 ^ 2 and D v (kFa,T/T F ) = (Aq) 2 (M v )qq . The normalization integral is 



^ = / dx/ (x)(l - f (x))x 6 = z I dx x 6 f + , (38) 



where z = e M//T is the fugacity and the collision matrix elements M^ t are given by 

_ f d 3 xd 3 y mTa 2 y ~ ~ 

Ms * " 7 l + 2mTaV 2 " 2 8 ' { ) 

with = f (xi)f (x 2 )(l - fo(x 3 ))(l - fofa)) and 

C t = [C t ( Xl )x\ 3 + C t {x 2 )4 - C t (x 3 )4 j - Ctixjxlt] /2 . (40) 

As in the case of thermal conductivity we can analytically compute D v by restricting the calculation 
to the lowest order orthogonal polynomial, and by considering the high temperature limit. We find 



f 5VmT 
J 32a 2 0? 


a 


| 15(mT) 3 / 2 
[ 320F 


a 



^ = s ( 41 ) 



In the superfluid phase the shear viscosity is dominated by binary phonon scattering. The leading 



low temperature behavior in the unitarity limit was calculated in 19]. Their result is 



3^7/2 (Tf 



T) H = 5.06 x 10- V 11 l-f ) Vo- (12) 



B. Prandtl Number 



The Prandtl number Pr = cpr]/(pK) characterizes the relative importance of thermal and mo- 
mentum diffusion. Here, p is the mass density of the fluid and cp is the specific heat at constant 
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pressure. In the high temperature limit the specific heat can be calculated using the Virial ex- 
pansion, and at low temperature it can be computed using the effective phonon lagrangian. The 
results are summarized in Appendix [B] The leading term at high temperature is determined by 
the equation of state of an ideal gas of spin up and down atoms. In this limit we get c p = 5n/2 and 
p = mn. The high temperature limits of k and 77 are given in Eq. (127)) and Eq. (j4ip . Putting these 
results together we get 

Pr = 2/3 T > Tp (43) 

independent of the size of the scattering length. This result is remarkable because it implies that 
the Prandtl ratio agrees with the classical gas result even though the shear viscosity and thermal 
conductivity are not classical. This result can be traced to the fact that in the large T limit the x 
and y integrals in Eq. ( l24"j) and (139)) factorize. The interaction only depends on y (related to the 
momentum transfer) and the leading y dependence of B\ ■ Bi and (% • Co is identical. This implies 
that both 77 and k can be expressed in terms of the same transport cross section, and that there is 
no dependence on the scattering length in the Prandtl ratio. In Fig. [5] we plot the Prandtl number 
as a function of T/Tp for negative values of kpa. As discussed in more detail in Appendix [B], the 
Virial expansion is smooth across the BCS/BEC transition, but our approximation for the bound 
state energy of the atom-atom dimer, Ep = l/(2ma 2 ), breaks down as we go to far towards the 
BEC side of the transition. For this reason we only show the Prandtl number for negative values 
of kpa. We observe that the Prandtl number increases significantly as the temperature is lowered, 
and that it is only weakly dependent on kpa. 



V. SOUND ATTENUATION 



Using the results for the shear viscosity and the thermal conductivity we can also compute the 
sound attenuation coefficient. The intensity of a plane wave decreases as e~ 2lx where 7 = aw 2 is 
the absorption coefficient and u is the frequency of the wave. In the normal phase the parameter 
a is given by [35] 



a 



1_ 

u 2 



2pc 3 



(i 


- L )] 


\c v 


Cp) 



(44) 



where p is the mass density of fluid, c is the speed of sound and ( is the bulk viscosity. At unitarity 
C = and we can write 



a 



2r] 
3pc 3 



1 + 



3 Ac P 

4~p7 



10 / * , * 



(45) 
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FIG. 5: The left panel shows the Prandtl number as a function of T jTp for different values of k F a. The 
right panel shows the Prandtl number as a function of kpa for several values of T/Tp. The dashed line 
shows the high temperature limit Pr = 2/3. 

where Acp = (c> — cy)/cy and rj = (mTp) 3 / 2 as in Eq. (l37|) . The quantities a* and a* are plotted 
in Fig. [6j 

In the high temperature limit the Prandtl number is Pr ~ 2/3 and Acp ~ 2/3. This implies 
that that the contribution to the sound attenuation coefficient due to thermal conductivity is 3/4 
of the contribution due to shear viscosity. Also, both Pr and Acp are only weakly dependent on 
the scattering length, see Fig. El As the temperature is decreased the sound absorption coefficient 
is increasingly dominated by shear viscosity. 

Below T c there are two sound modes, ordinary (first) sound and second sound. The absorption 
coefficient for first sound is analogous to Eq. f|4*4l) . except that the contribution from k is absent. 
This implies that damping of first sound is entirely due to shear viscosity. The absorption coefficient 



of second sound is given by [36l 



2 pc^ p n [3 p s C P 



(46) 



where p nyS are the normal and superffuid densities, C2 is the speed of second sound and Q are the 
four different coefficients of bulk viscosity in the superffuid. Onsager's symmetry principle implies 
C4 = Ci> an d conformal symmetry requires that at unitarity Ci = C2 = 0. The remaining coefficient, 



(3, vanishes as £3 ~ T 3 as the temperature goes to zero 21J. Using 77 ~ 1/T 5 and k ~ T 2 we 



conclude that the damping of second sound is also dominated by shear viscosity. 
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FIG. 6: a* and a* as a function of T/Tp. Note that the ratio of the two approaches a constant (3/4) as 
we go to larger values of T jTp. 

We should note that at low temperature the phonon mean free path becomes very long, and 
the normal fluid ceases to behave hydrodynamically. A rough estimate of the relaxation times for 
shear viscosity and thermal conductivity is t v ~ i]/(nT) and r K ~ mn/(cpT). This implies that 
the relaxation time for shear viscosity grows very quickly, r rj ~ 1/T 6 , whereas the relaxation time 
for thermal conductivity grows more slowly, r K ~ 1/T 2 . 

VI. CONCLUSIONS 

In this work we studied the thermal conductivity of a dilute Fermi gas. Our main results are 
summarized in Fig. (jl]). We observe that the thermal conductivity scales as n ~ m 

i/2 T s/2 in the 

high temperature limit and as k ~ m^T^/ 2 (T jTp) 2 in the low temperature limit. The Prandtl 
ratio equals 2/3 in the high temperature limit, and increases as the temperature is lowered. This, 
together with the fact that Acp = (cp — cy)/cy decreases at low temperature, implies that the 
sound attenuation length in the strongly coupled normal fluid regime T c < T < Tp is dominated 
by shear viscosity. In order to assess the feasibility of measuring the sound attenuation length we 
consider the ratio of wavelength A to the attenuation length 7 -1 . We find 



\ = 4?r ( k \(V 
7 v/SfUjJU 



3 Ac P 
1 + 4~PT 
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k F J \n 



where we have used £ ~ 0.44 and Acp/ Pr ~ 0. A rough estimate of the minimum value of 77/n can 
be obtained by extrapolating the kinetic theory result down to T = T c . We get 

rj\ 45vr 3 / 2 / T, \ 3/2 



0.15. (48) 



.nJmin 64a/2 V^V, 

Finally, we can express k/kp in terms of the wavelength and the trap parameters. We obtain 

k R z vr/3 2 / 3 



(49) 



k F A (3A0 1 / 3 ' 

where R z is the longitudinal cloud size, (3 = R±/R Z is the trap asymmetry, and A" is the number 



of particles. For the experiment reported in 



22| we get k/k F ~ 0.046 and A7 ~ 0.08. This implies 



that a measurement of the sound attenuation length is within reach. Higher sensitivity can be 
achieved by using more elongated traps or smaller wavelengths. 

There are a variety of improvements that need to be made once experimental data become 
available. The most important is to replace the plane wave approximation which is the basis of 
Eq. (j47p by a more detailed calculation of the wave profile in a given trap geometry. It will also be 
interesting to study attenuation of first and second sound in the superfluid phase. 
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discussions. 



Appendix A: Trial Functions 

In this appendix we provide explicit expressions for some of the orthogonal polynomials and 
associated normalization constants. The orthogonality condition for the polynomials B^ F,H \x 2 ) is 
(see Eq. (HQ) 

Ai F ' H) S st = J dxf (x)(l ± h{x))x'B^ H \x 2 )B^ H \x 2 ) , (Al) 
where Bq F ' H \x 2 ) = 1 and 

I^V)=^ + X>if^ (A2) 

i=0 

(F H) 

for s > 1. We will determine the coefficients a si ' iteratively. We define 

= J dx /„(*)(! ± f (x))x 4+n , (A3) 



16 



so that A = Co- We also find 



a w = -c 2 /co, (A4) 
a 2o — ( c 2Ce — cl) / (cqC^ — cl) , 

a 21 — (CoCq — C2C4)/ (C0C4 — Cg) , 

and A\ = C4 — Cg/co- The integrals c n depend on the quasi-particle dispersion relation. In the case 
of the phonons with E p = vp, we can evaluate A\ to be 

In the case of the fermion with E p = p 2 /2m and a chemical potential, we can only evaluate A\ in 
the case of large temperature which gives 

A{ * z 1 ^- (A6) 

where z = e^ T is the fugacity. 

The orthogonal condition of the shear viscosity's trial function polynomials g(x) = Yl c sC s (x 2 ) 
is different from the one of thermal conductivity. Here, the polynomial is starting from s = term 
since this set of polynomials satisfies the constraint of momentum conservation law automatically 
(due to the tensor nature of the kernel for shear viscosity). The orthogonality condition is then 

A v s 5 st = j dx/ (x)(l - f (x))x 6 C s (x 2 )C t (x 2 ) (A7) 

where the trial function C s (x 2 ) = x 2s + YliZo a siX 21 . By the orthogonality condition Eq. ( 1A7j) and 
the same recursion method that we used before, it is easy to get the coefficients a S i and also A$ 
which is given by 



,.2 



Al = j dz/o(s)(l - h{x))x & = z j dx f + * z 1 ^- (A8) 



where the analytic result is valid in the high temperature limit. Note that Af = A^ in the high 
temperature limit; however, this is a mathematical coincidence coming from the similarity of the 
different integrals and the properties of the T-function. 



17 



Appendix B: Thermodynamic Properties 



1. Virial Expansion 

In the high temperature limit, T ^> T F , the pressure can be expanded in powers of the fugacity 
z = e M / T . This expansion is known as the Virial expansion. We write 

P{^T) = U —'£b l z\ (Bl) 
1=1 

where A = a/27t /mT is thermal wave length an v is the spin degeneracy factor. The first Virial 
coefficient isjh = 1. The second Virial coefficient was first calculated by Beth and Uhlenbeck 38], 



see also 



39 



40| | . Here we use these results in order to compute Acp. There are two contributions 
to b 2 . The first is due to quantum statistics, which becomes important as n\ 3 ~ 1, the other is 
related to the atomic interaction. The correction due to quantum statistics is 

Bosons , 




(B2) 

Fermions . 

The interacting contribution to the second Virial coefficient is 

b 2 -b° 2 = V2j2(e~ tn/T ~ e" e ° /T ) , (B3) 

n 

where e n is the energy of the n'th two-particle state in the center of mass frame and e° is the 
corresponding non-interacting energy. In the universal regime these energies are determined by the 
scattering length only. We get [40] 



b 2 -bl 



-^e x [1 - erf(|x|)] x < 
V2e x2 - -^e x2 [1 - erf (x)} x>0 



where erf(x) is the error function, x = \/(a\/2rt) and the bound state energy in the x > case was 
assumed to be Eb = I /(ma 2 ). We can write down the asymptotic forms in the limits of zero and 
infinite scattering lengths as 



b 2 -b° 2 



^2 e l/{rnTa*) a ^ + , (B4) 

-4(l + ^rr) a^±oo. 

K V2 K n \a\' 
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FIG. 7: The left panel shows cy in units of uq as a function of T/Tp no- Here, uq = uz/X 3 is the density 
of the free gas. The right panel shows the temperature dependence of cp. 

The interaction part of 62 approaches a finite value as a — > 00, and is of the same order as the 
effects of quantum statistics. One can also see that the total 62 will start off negative for small and 
negative scattering length, but eventually becomes positive as a — » — 00. This implies that there is 
a value of the scattering length for which the effects of the interaction and of quantum statistics 
cancel. For positive scattering lengths, the contribution of the interaction is positive and always 
larger than the quantum correction, so that 62 > for all a > 0. In the BEC limit (a — > + ) 
the binding energy of the bound state is large and 62 becomes large also. In this limit, however, 
universality breaks down and the binding energy cannot be expressed in terms of the scattering 
length only. 

2. Specific Heat 

The specific heat can be expressed in terms of partial derivatives of the pressure. The inde- 
pendent thermodynamic variables are (/i,T). The first derivatives with respect to theses variable 
determine the entropy density and pressure 



dP 



dP 



(B5) 




dfj, T 
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In order to compute the specific heat at constant volume we use V = N/n and write 



Cy = T 



ds 
df 



d(s,V) _ d(s,V)/d(T^) 
v d(T,V) d(T, V)/d(T, /i) 



T 



ds 
df 



[{dn/dT)U> 
(dn/djj)\ T 



(B6) 



where the Jacobian is defined as 



d(a,b) I d a /d c \d 9a/dd\ 



d(c, d) 



db/dc\. db/dd\ 



In order to compute cp we make use of the relation between cp — cy and the thermal expansion 
coefficient (l/V)(dV/dT)\p. This relation is given by 



cp — Cy = — T 



\{dV/dT)\ P f 
(dV/dP)\ T 



The partial derivatives are 



dV 




s dn 




dn 




dV 




dn 


dT 


p 


n dfi 


T 


dT 


1-1 


dP 


T 





which gives 



cp = Cy + T 



%(dn/dfj) 


T - (dn/dT) 


- 2 


(dn/dfi) 


T 



(B7) 



(B8) 



(B9) 



a. High Temperature Phase 



In the high temperature phase we use the Virial expansion of the pressure including terms up 
to second order, P = ^(z + b 2 z 2 ). Using the results of the previous section we find 



vz 

CV= V 



3 15 db 2 2 d 2 b 2 

2 + ^ Zh ' zT df + zT d^f 



(BIO) 



cp 



vz 

cv+ p 



1 + 5zb 2 - 2zT^ 



We recover the result for a classical non- interacting gas by setting b 2 = and n = p. We show 
the result for an interacting gas in Fig. [71 We observe that the results are smooth across the 
BCS/BEC transition. We note, however, that on the BEC side the Virial expansion breaks down 
at temperatures larger than the Fermi temperature. This is related to the presence of a deeply 
bound state. 
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b. Superfluid Phase 



At zero temperature the pressure is given by 



^ = I^™ 3 V' 2 , (bid 



and the density is 



_ dP 



j, ~ -^w" ■ (B12) 



The leading low-temperature correction arises from phonons. The phonon contribution to the 
pressure is 

vr 2 T 4 /3m\3/2 



ft = -9o-Ur.» (B13) 



Using these results we find 

2tt 2 T 3 /3m\3/2 



f^V' Z + 0(T'), (B14) 



15 V2/i 
c P = c y + £>(T 7 ) . 

These results imply that Ac p vanishes at low temperature much faster than cp and cy. 
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